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Abstract -We study the optimal routing on multilayered communication networks, which are 
composed of two layers of subnetworks. One is a wireless network, and the other is a wired network. 
We develop a simple recurrent algorithm to find an optimal routing on this kind of multilayered 
network, where the single-channel transmission mode and the multichannel transmission mode 
used on the wireless subnetwork are considered, respectively. Compared with the performance 
of the shortest path algorithm, our algorithm can significantly enhance the transport capacity. 
We show that our methods proposed in this letter could take advantage of the coupling of the 
two layers to the most extent, so that the wireless subnetwork could sufficiently utilize the wired 
subnetwork for transportation. 



I. Introduction. — Transportation is an important 
problem in natural and engineering systems. The study 
of transportation optimization on complex networks has 
attracted a lot of interest in the past decade [Tl[2] . These 
studies include the influence of network structure on the 
transportation [SHZl, features of different transportation 
models [8l-[T3] , and offer strategies to improve transporta- 
tion from different viewpoints [HIITS]. The shortest path 
route, which takes the shortest-distance paths between the 
origins and the destinations, is a most common strategy 
for routing. However, this approach easily leads to con- 
gestion on some popular nodes having the most number 
of routes passing by. Therefore, to overcome this short- 
coming several efficient strategies have been developed 
in a static scheme by focusing on the rerouting of con- 
gested flow so as to lead the system from congestion to 
free flowing [T6l[T7|. It is proved that to find the optimal 
routing of a transportation network, which has maximum 
transportation capacity without congestion, is an iVP-hard 
problem. It is remarkable that Bassler et al [TTj recently 
presented routing strategies for wired networks and wire- 
less networks respectively for the purpose of increasing 
the transportation capacity of these networks, which out- 
perform many other methods. Their algorithm is a sub- 
optimal solution and takes the running time 0{N^ logN). 
Hence, the transportation capacity can be much improved 



if one can judiciously choose the routes. 

With the rapid development of communication net- 
works, wired and wireless networks have been widely ap- 
plied in many aspects of society, which include the Inter- 
net, which is based on the wired infrastructure and mobile 
telecommunication, which is based on the wireless frame- 
work. Despite the works introduced above, relatively little 
attention has been paid to the integration of the wired and 
wireless networks. The issue of taking the wired and wire- 
less networks as a whole subject has become increasingly 
fundamental due to the widespread overlap of their appli- 
cations. In this paper, we introduce an intuitive model 
which captures the important characteristics of the wired 
and wireless networks, and then we study transportation 
optimization on such integrated networks. Specifically, we 
first construct a wired network and a wireless network. 
Then, we choose the same number of nodes in each net- 
work and combine them pairwise to form a new one. Thus, 
the previous wired and wireless networks become the two 
layers of the new network. Notice that some characteris- 
tics of the wireless network are different from the wired 
one. One is the absence of the effects of rich hub connec- 
tivity and small shortest path length, which by contrast 
are important features of wired networks. The other is a 
constraint in the process of packet transmission. The con- 
straint comes from the broadcast transmission mode of the 
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wireless device. When a wireless device is sending packets, 
other devices within the one- hop transmission range of the 
sender or the receiver are forbidden to send or receive in- 
formation packets if they arc working in the same channel, 
which is called medium access control. However, the wired 
networks do not have such constraint. After that, wc pro- 
pose a heuristic routing method to implement the efficient 
routine on the multilayered networks by considering the 
different features of the wired and wireless networks |17j . 
Comparing with the shortest path method, wc find our 
method can significantly enhance the transportation ca- 
pacity of the two-layered network. 

The remainder of this paper is structured as follows. In 
sec. II, we present the model of the multilayered commu- 
nication network. In Sec. Ill, we present our algorithm 
for single channel transmission mode and show the results 
of the algorithm. In Sec. IV, we study the general situa- 
tion where the wireless networks work in the multichannel 
mode. In Sec. V, wc give our conclusion. 

II. The Model. — We initially construct two sub- 
networks. One. called wired layer, is a wired network with 
iVw nodes and the other, called wireless layer, is a wireless 
network with iVL nodes. The wired layer is constructed 
with the Positive-Feedback Preference (PFP) Model [18] 
which accurately reproduces many topological properties 
of AS-level Internet. The wired layer constructed by the 
PFP model has the effects of rich hub connectivity and 
small shortest path length, and has an average degree 
~ 5.4. For the wireless layer, we use the minimum de- 
gree geometric network model |19| . which is appropriate 
in describing wireless networks. In the random geomet- 
ric network model, a pair of nodes can make contact with 
each other when their distance is shorter than a critical 
value. However, this critical value needs to be provided 
by an external control authority which does not exist in 
ad hoc networks. For the minimum degree geometric net- 
work model, all the nodes are randomly distributed on a 
square and they can decide on their own contact radius 
by increasing their contact power until they have at least 
fcmin mutual neighbors. Therefore, with a given kmin, all 
the nodal degrees in the wireless layer will be no less than 
fcmin. After generating the two layers, we randomly choose 
a node in the wired layer and a node in the wireless layer 
and then merge them into a single one. The newly merged 
node preserves all the connections that the previous two 
nodes have. Thus, this node could serve as an interface of 
the two layers. We continually perform this process until 
A^i pairs of nodes from the two layers are merged to form 
Ni nodes. For simplicity, we refer to these Nj nodes as 
interfacing node which belong to the both layers, and we 
refer to the nodes in the wired (wireless) layer other than 
the interfacing nodes as wired (wireless) node. Thus, A^w 
(A'l) is the number of nodes in the wired (wireless) layer 
and the size of the whole network is A^ = N-w + A'l — A'l. 
Figure [1] offers an illustration of the network structure. In 
this paper, we set the minimum degree fcmin as 8. We ex- 
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Fig. 1: Schematic illustration of the multilayered communi- 
cation network with broadcast constraint, where node A, B 
and C are wireless nodes, node D is an interfacing node, 
and node E is an wired node. The dashed lines denote 
the wireless connections and the straight lines denote the 
wired connections. In this figure, the effective betweenness-to- 
capacity ratio of node A, denoted as {B/CYa , is (B/CYa = 
{B/C)a + {B/C)b + {B/C)c + (B/C)d, and correspondingly 
(B/C)"/ = {B/C)d + {B/C)a and (B/C)f = (B/C)e. 



amined other values of fcmin and find that the results are 
robust. In the simulations, all the averaged results and 
their standard deviation (which is indicated by the error 
bars), are obtained from 500 different realizations, if not 
otherwise specified. 

III. Single Channel Mode. — Wireless networks 
have a broadcasting constraint when the nodes work on 
the single channel mode [TH] • The constraint is that when 
a node is sending an information packet, all its incom- 
ing links are forbidden to send packets in order to avoid 
packet collisions. This constraint causes a big obstruction 
for the efficient information transmission. In the transmis- 
sion process, the wireless nodes work on the first-in-first- 
possible-out basis. That is, a wireless node attempts to 
send its first-in-line packet. When the intended recipient 
is blocked due to the broadcasting constraint, the node 
then tries to send its second-in-line packet and so on, un- 
til cither an idle recipient is found or the end of the queue 
is reached. While, for the wired networks there is no such 
constraint, and therefore the wired nodes simply work on 
the first-in-first-out basis. The interfacing nodes adopt the 
transmission process according to whether the intended 
recipient or the sender is a wired node or a wireless node. 

When a new packet is added to the network at a node, 
it is appended at the end of the queue of this node. A 
packet is removed when it reaches its destination. When 
the routes between all pair of nodes are given, the be- 
tweenness of a node i is defined as Bi = Yis^tPi(st) /p(st), 
where p(st) is the number of routes going from node s to 
node t and Pi(st) is the number of routes going from node 
s to node t and passing through node i. Wc assume that 
a node i has the processing capacity of handling Ci in- 
formation packets per time step and different nodes may 
have different processing capacities. Thus, the value of 
the betweenness-to-capacity ratio of node i, denoted as 
(B/C)i, equals to Bi/d. 

An important consequence of the broadcast constraint 
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is that a node in the wireless layer needs not only to pro- 
cess the information packets on itself but also to process 
the packets on its wireless neighbors. This situation leads 
to a modified version of the bctweenness-to-capacity ra- 
tio [H], which is the summation of the betweenness-to- 
capacity ratios of the objective node and its wireless neigh- 
bors. However, the wired nodes do not need such modifi- 
cation. Hence, the effective betweenness-to-capacity ratio, 
denoted as (i?/C)®*, for different nodes is different. The 
effective betweenness-to-capacity ratio of a wireless node is 
equal to the sum of its own betweenness-to-capacity ratio 
and the betweennesses-to-capacity ratio of all its neigh- 
bors. The effective betweenness-to-capacity ratio of an in- 
terfacing node is equal to the sum of its own betweenness- 
to-capacity ratio and the betweennesses-to-capacity ratio 
of its neighbors connected with its wireless connections. 
Finally, the effective betweenness-to-capacity ratio of a 
wired node is equal to its own betweenness-to-capacity ra- 
tio (see Fig.[T]). The maximum transportation capacity de- 
pends on the maximum effective betweenness-to-capacity 
ratio which is denoted as (B/C)^^^. The optimal routing 
is achieved when {B/C)'^^^ is minimized. Adopting the 
algorithm that is used for single layered network [T7], we 
extend it to the multilayered communication network to 
reduce the (-B/C)^^ through the following steps: 

1. Assign every link a unit weight and compute the 
shortest path between all pairs of nodes and the between- 
ness of every node. 

2. Calculate the effective betweenness-to-capacity ratio 
of all the nodes and find the node which has the highest 
score, (a) If the node is on the wireless layer, increase 
the weights of all the links sourced from this node and 
its neighbors by half unit weight, (b) If the node is an 
interfacing node, increase the weights of the links sourced 
from this node and its neighbors in the wireless layer by 
half unit weight, (c) If the node is on the wired layer, 
increase the weight of the links sourced from this node by 
half unit weight. 

3. Recompute the shortest paths and the bctwccncsscs. 
Go back to step 2. 

4. The algorithm stops when the value of the maximum 
effective betweenness-to-capacity ratio is stable. 

We note that the weight of a link can be regarded as the 
measure of the length of this link. Hence, a larger weight 
corresponds to a longer distance, and a packet will less 
likely pass through the corresponding link because of the 
higher cost required. 

A plot of the maximum effective betweenness-to- 
capacity ratio {B/C)'^^^ obtained from the shortest path 
method and the optimal algorithm versus network size is 
shown in Fig. [2ja). Obviously, (B/C)^^ increases with 
the increasing of the size of the networks. The inset shows 
the ratio p of the results obtained from the optimal algo- 
rithm to those obtained from the shortest path algorithm. 
We observe that when the network size is small, the effect 
of the optimal algorithm is very limited. This is because 
of the fact that those small networks do not leave much 




Fig. 2: (a) The maximum effective betweenness-to-capacity ra- 
tio for the shortest path algorithm (solid symbols) and for the 
optimal algorithm (open symbols) as a function of A^v and A^'l. 
In each network, Ni — 10 and A'w = Nh. The inset shows the 
ratio p of the results obtained from the optimal algorithm to 
those obtained from the shortest algorithm, (b) The value of 
(_B/C)^^ after performing the optimal algorithm as a function 
of the number of interfacing nodes for different network sizes. 
The inset is a log-log plot of the main plot. In this figure, all 
the nodes have the capacity C = 1. 

space for the algorithm to adjust the transmission routes 
to reduce the (B / C)'^^^. However, with the growth of the 
network size, there could be more alternative routes for a 
transmission. Therefore, it could be easier for the algo- 
rithm to find a route without passing the node that has 
the maximum betweenness-to-capacity ratio, which results 
in a considerable decrease in the (iJ/C)^^;^^- Hence, this 
algorithm could be more effective for a larger network, 
which is a good point as real communication networks are 
generally very large. 

The interfacing nodes take on an important role for the 
whole network. Since the routes of the inter-layer com- 
munication need to pass through the interfacing nodes, 
when the number of interfacing nodes is small they may be 
the bottleneck. Figure [Ub) shows the maximum effective 
bctwccnness-to-capacity ratio {B/C)°^^^ from the optimal 
algorithm as a function of A'l for different network sizes. 
We can see that when Ni is small the value of [B/C)'^^^^ 
can be very large, which means the algorithm has limited 
effect in this situation because of the bottleneck effect of 
the interfacing nodes. However, with the increase of the 
Ni the value reduces quickly and the situation gets much 
improved. 

IV. Multichannel Mode. — Single channel mode is 
the primary mode for wireless communication, but its effi- 
ciency is low. Thus, this mode is not suitable for large net- 
works. In fact, the multichannel mode is widely adopted 
in engineering, and explicit regulations have been devel- 
oped for proper operation |20| . A common way is to let 
wireless devices consistently change their working channel 
according to a pre-encoded random sequence. In this way, 
when a group of nodes located nearby are transmitting in- 
formation, the possibility of them using the same channel 
is low. Even when this happens they can simply turn to 
other channels to avoid such situation. Thus, the situation 
in which different pairs of communicating nodes happen to 
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Fig. 3: (a) (i3/C)max as a function of A^i for the shortest path 
algorithm (black square symbols), the optimal algorithm for 
Random Placement (RP) (solid circle symbols), the optimal 
algorithm for Optimal Placement (OP) (open circle symbols) 
and the theoretical lower bound which is {2Nw Ni,) / (NiC) . (b) 
Average transmission distance between wireless nodes cJll as 
a function of Ni for the shortest path algorithm (solid square 
symbols), the optimal algorithm for RP with C = 1 (solid 
circle symbols), the optimal algorithm for RP with interfacing 
nodes having C = 4 (solid triangle symbols), and the optimal 
algorithm for OP with interfacing nodes having C — 4 (open 
triangle symbols). The inset shows the (-B/C)max of the cases 
that the interfacing nodes have C = 1 and C = 4 as a reference. 
All the networks in this figure have Aw = 200 and Al = 1000. 
All the nodes other than the interfacing nodes have C = 1. 

work on the same channel is well avoided. Since nodes that 
work on different channels do not influence each other, 
the effect of the broadcast constraint can be ignored in 
the multichannel mode and the communication efficiency 
could be significantly enhanced. In the following, we study 
the transport optimization of the multilayered network in 
which the wireless layer works in the multichannel mode. 

Since the broadcasting constraint can be ignored under 
the multichannel mode, it is not necessary to use the mod- 
ified version of betweenness-to-capacity ratio as did in the 
previous section. In this mode, the effective betweenness- 
to-capacity ratio is just the betweenness-to-capacity ratio 
of the node itself. Thus, in the multichannel mode the 
algorithm works as follows: 

1. Assign every link a unit weight and compute the 
shortest path between all pairs of nodes and the between- 
ness of every node. 

2. Calculate the betweenness-to-capacity ratio of all 
the nodes and find the node which has the highest value. 
Increase the weight of the links reached by this node by 
half unit weight. 

3. Recompute the shortest paths and the betweennesses. 
Go back to step 2. 

4. The algorithm stops when the value of the maximum 
betweenness-to-capacity ratio is stable. 

An interesting point of the multilayered network is 
the influence of the interfacing nodes. Since all the 
inter-layered transmission need to go through the inter- 
facing nodes, we may estimate the lower bound of the 
maximum betweenness of the interfacing nodes. From 
the definition of the betweenness, we have ~ 



SiGi 12s=it Pi{^'^) I P{^^) ^ where I is the set of the interfac- 
ing nodes. This equation could be further separated into 
four cases. Those are (i) s G W, t £ W; (ii) s £ C, t € £; 
(iii) s e yy, t G £; and (iv) s e C, t £ W; where W and C 
are the sets of the wired nodes and the wireless nodes, re- 
spectively. Since each packet for inter-layered communica- 
tion passes through the interfacing nodes at least once, for 
cases (iii) and (iv) we have ^ligi J2sgw tecP^i^'^)/p(^'^) ^ 
NwNi^ and E^ei Ese£,tewK(s*)MsO > ^wA^l- The 
equal sign is obtained when each packet for inter-layered 
communication passes through the interfacing nodes only 
once. Therefore, even all the routes for intra-layered com- 
munication are restricted in the respective layer, that is 
Pi£x{st) ^ for cases (i) and (ii), we have J2iex — 
2AwA^L- Since there are A^i interfacing nodes, the max- 
imum betweenness of the interfacing nodes should be no 
less than 2Nt^Ni^/Ni. This value is reached when the in- 
terfacing nodes take the same amount of the information 
packets transmitted between the layers and routes of the 
intra-layered communication are confined in that layer. 
Thus, when A^i is small, the largest betweenness of the 
interfacing nodes can be very large. Figure ^{a) shows 
the maximum betweenness-to-capacity ratio {B/C)ma.x as 
a function of A^i obtained from the shortest path algorithm 
(square symbols), the optimal algorithm (circle symbols), 
and the theoretical lower bound of the interfacing nodes 
(2AfwAfL)/(AfiC) (black curve), for the case of C = 1, 
A^w = 200 and A^l = 1000. Similarly to the single chan- 
nel mode, (i?/C)max decreases with the increase of Ni, 
which means there also exists the bottleneck effect in the 
multichannel mode when the number of interfacing nodes 
is small. Moreover, the circle symbols almost overlap with 
the black curve when A^i < 6 and are just a little larger 
when A'l > 6, which shows that the optimal algorithm 
can reduce the {B/C)max nearly to the theoretical lower 
bound and thus to the optimal level. We note that when 
A^i = 1 the values of (-B/C)max of the two algorithms are 
very close. This is because when there is only one inter- 
facing node, this node has to handle all the inter-layer 
communications whether or not the optimal algorithm is 
adopted, which causes the node to have the (i?/C)niax and 
the value is approximated to {2NwNi,)/{NiC) = 4 x 10^. 
Furthermore, for the optimal algorithm when A^i = 2 the 
value of (i?/C) max drops by around a half compared to 
the case of A^i = 1, which means this algorithm can well 
distribute the information packets to the two interfacing 
nodes. Furthermore, we observe that when A^i is small the 
distribution of {B/C)ma.x obtained from the optimal algo- 
rithm under different realizations is squeezed in a narrow 
range. For example, for A'l = 1, the ratio of the standard 
deviation of the distribution of (-B/C)niax with performing 
the algorithm to that without performing the algorithm 
is smaller than 0.01. This phenomenon is distinct from 
what is observed in the single channel mode. In the sin- 
gle channel mode, the distributions of the (B/C)^^^ with 
and without performing the optimal algorithm spread in 
the ranges with comparable sizes, and both can be well 
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fitted by the Gumbel distribution [3T] . The reason of this 
phenomenon is that in the multichannel mode, the opti- 
mal algorithm could always tune the (-B/C)inax close to 
the theoretical lower bound resulting in the fact that the 
values of (i?/C)niax for different realizations are very near 
to each other, which further proves the effectiveness of the 
algorithm. 

As the maximum betweenness-to-capacity ratio of the 
interfacing nodes should be no less than (2A'w-^l) / (NiC), 
when the value of (i?/C)max is close to this lower bound, 
the routes for the intra-layered communication need to be 
mainly restricted in the respective layer to avoid going 
through the interfacing nodes. Thus, two wireless nodes 
cannot use the wired layer as well as its effects of rich 
hub connectivity and small shortest path length to shorten 
their distance in this situation. In order to improve the 
transmission efficiency by reducing the transmission dis- 
tance, without increasing the value of (i?/C)inax, an effec- 
tive method is to increase the processing capacity of the 
interfacing nodes. To illustrate the effectiveness of this 
method we show the average transmission distance among 
all pairs of wireless nodes, denoted as ^ll, in Fig.|3](b). In 
this figure, d^L obtained from the shortest path algorithm 
as shown by the solid square symbols is the smallest one 
among all the cases, while cIll obtained from the optimal 
algorithm where all the nodes have C = 1 as shown by the 
solid circle symbols is much greater especially for large Ni. 
However, when we increase the processing capacity of the 
interfacing nodes to C = 4 while other nodes still have 
C = 1, the dLL has an obvious drop as shown by the solid 
triangle symbols which shows that increasing the process- 
ing capacity of the interfacing nodes could be of benefit to 
the wireless nodes utilizing the wired layer to reduce the 
transmission distance. 

Moreover, we consider that if the positions of the in- 
terfacing nodes in the wireless layer are well-distributed, 
the wireless nodes could access their nearest interfacing 
nodes with the least of hops on average. Thus, this place- 
ment could further shorten the transmission distance when 
wireless nodes utilize the wired layer for transmission. For 
simplicity, we call such placement of the interfacing nodes 
in the wireless layer as Optimal Placement (denoted as 
"OP"), and correspondingly the placement in which the 
interfacing nodes are randomly distributed in the wireless 
layer as we did in the previous part of the work is called 
Random Placement (denoted as "RP" ) . To obtain the OP, 
we first get A^i optimal positions with the algorithm intro- 
duced in the appendix. Then, in the stage of constructing 
multilayered network we choose Ni wireless nodes who are 
nearest to the A^i optimal positions as the wireless candi- 
dates to compose the Ni interfacing nodes. In this way, 
the composed interfacing nodes have OP in the wireless 
layer. The open triangle symbols in Fig. |3](b) show the 
results of the OP where the interfacing nodes have C = 4. 
We observe that the d^L is further reduced in this case. 
The inset shows the corresponding (i?/C)max for the cases 
that the interfacing nodes have C = 1 and C = 4. Obvi- 



ously, the (i?/C) max for the case of C = 4 is always smaller 
than that of C = 1. We note that (i3/C)niax of OP and 
RP are very close when the capacity of the interfacing 
nodes in both cases is the same. An example is shown 
in Fig. |3Ka) where the open circle symbols correspond to 
OP and solid circle symbols correspond to RP. Therefore, 
by combining the methods of using the optimal algorithm, 
increasing the capacity of the interfacing nodes which are 
a small fraction of the whole network, and the OP, both 
the transportation capacity and average transmission dis- 
tance can be much improved. Since these two factors are 
generally conflicting with each other, it is significant that 
both of them can be enhanced at the same time. 

V. Conclusion. — In summary, we propose a model 
to describe multilayered communication networks. In the 
model, there are two layers. One, called wired layer, is 
composed of nodes connected with wired links, and the 
other, called wireless layer, is composed of nodes con- 
nected with wireless links. The layers are connected by 
merging some pairs of nodes in different layers into single 
ones, called interfacing nodes. Then, we present a recur- 
rent algorithm to find optimal routing for this multilay- 
ered network for two different cases in which the wireless 
nodes can work in single channel mode and multichannel 
mode, respectively. In this algorithm, by gradually in- 
creasing the weight of the connections pointing to the node 
which has the maximum effective betweenness-to-capacity 
ratio, the burden on this node could be gradually reduced. 
By repeating this process, the final maximum effective 
betweenness-to-capacity ratio could be much suppressed 
and the transportation capacity could be enhanced signif- 
icantly. Our results show that this algorithm is an effec- 
tive navigation technique for multilayered communication 
networks. Since for the network on one hand the wireless 
layer may utilize the wired layer for transportation, on the 
other hand the interfacing nodes can easily become bot- 
tleneck, how to utilize the advantage of the multilayered 
networks while avoiding the restriction of the interfacing 
nodes is crucial for enhancing the transportation capacity. 
Our routing method proposed for the multilayered net- 
works can squeeze the interfacing nodes to the extent that 
their effective betweenness-to-capacity ratio is almost the 
same, so that the bottleneck effect is suppressed to the 
most extent. By further cooperating with the methods 
of increasing the processing capacity and the OP of the 
interfacing nodes, the wireless layer may have a sufficient 
utilization of the wired layer for transportation. Hence, 
our method could serve as a good candidate of the stan- 
dard technique for optimal transportation on multilayered 
communication networks. 

Appendix. A method of placing an arbitrary 
number of nodes on a square in a v^rell-distributed 
pattern. — Suppose there are N nodes to be dealt with, 
we first regard the N nodes as N solid circles with the same 
radius which are required to be confined in the square. It 
can be proved that placing the nodes on the square so 
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Fig. 4: (a) An example of A'' = 2, where labels 1 and 2 indicate 
the centers of circle 1 and 2. We denote the candidate radius 
of circle 1 and 2 as n and r-2, respectively. For circle 1, xl, 
a^R, 2/u and j/d are the distances between its center to the four 
boundaries of the square, and di2 is the distance between the 
two centers. Hence, r\ = minjiL, 2;r, yu, J/d, cii2/2}. Similarly, 
we can get the value of r-^. Thus, the value of rmin is rmin = 
min{ri,r2}. Dashed circles are those plotted by the candidate 
radii, (b), (c) and (d) show the examples of the cases A'' = 4, 
7 and 15, respectively. 

that they are well-distributed is equivalent to place these 
soUd circles so that they can have the largest identical 
radius, and the positions of the centers of these circles are 
the positions in vifhich to place the nodes. Following we 
present the algorithm of placing these solid circles so that 
they can have the largest identical radius: 

1. Initially, randomly distribute these nodes on the 
square. Then, assign each node a candidate radius which 
is the shortest one among the distances that are from the 
node to the four boundaries of the square or as half as 
those from the node to the other nodes. The identical 
radius of these solid circles which satisfies the above con- 
ditions, denoted as rmin, is the smallest one among all the 
candidate radii (see Fig. |3^a)). 

2. Disturb the position of a randomly chosen node and 
calculate the rmin again. If r-min decreases, the disturbance 
is discarded and the position of the disturbed node is re- 
covered. If Tmin is unchanged, the disturbance is accepted 
and the position of the disturbed node is updated. If rmin 
increases, the disturbance is accepted and the position of 
the node as well as the value of the r,nin are updated ac- 
cordingly. 

3. Go back to step 2 and perform a new disturbance 
until rmin is stable. 

After performing the algorithm, the final positions of 
the centers are the positions in which to place the nodes. 

* * * 
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